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Abstract 

Multijet events at large transverse energy (J2, Et > 420 GeV) and large 
multijet invariant mass (mjets > 600 GeV) have been studied by the CDF 
Collaboration at the Fermilab Tevatron. The observed jet multiplicity dis- 
tribution can be understood in a QCD inspired exponentiation model, in 
regions of phase space which require going beyond fixed order perturbation 
theory. 
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With the completion of run I at the Fermilab Tevatron, a sufficient amount of data has been 
collected to study rare hard scattering events, with jet transverse energies in the hundreds of 
GeV range. Such high Et jet events are interesting in their own right, as the discussion of a 
possible "excess" in the single jet Et spectrum above Et ~ 200 GeV [0] has shown. Another 
aspect is the fact that events with jets in this Et range, while rare at the Tevatron, will be 
produced abundantly at the much higher energy of the CERN LHC, where they will constitute 
important backgrounds to searches for new heavy degrees of freedom. Thus, the study of high 
Et QCD events at the Tevatron allows one to develop strategies for new physics searches at the 
LHC. 

The aspect to be considered below is the multiplicity of additional soft jets which arise via 
the emission of soft gluons in a hard scattering process. At the LHC, one application is the 
study of weak boson scattering. In events like qQ — > qQWW, which are mediated by t-channel 
exchange of electroweak bosons, medium Et jets in the central region are a rare occurrence. 
This is a consequence of color coherence between initial and final state radiation which leads 
to gluon emission mainly between the forward scattered quarks and the beam directions 0. 
At the same time the modest transverse momentum of the scattered quarks severely limits the 
transverse momentum of emitted gluons |^. Typical background events like ti — hW^hW~ or 
qq — > W~^W~ , on the other hand, show a large probability for QCD radiation in the central 
region, with gluon transverse momenta sufficient to produce visible jets 0]. As a result, a central 
jet veto is quite effective for background suppression 

In applications like these it is necessary to correctly model the angular distribution of emitted 
minijets and to reliably determine the hard scales which govern their transverse momenta. These 
features are automatically included by using full tree level QCD matrix elements. For an effective 
jet veto, however, one is interested in the phase space region where the probability of extra QCD 
radiation becomes of order unity, and this is exactly the region where fixed order perturbation 
theory ceases to be applicable. In Refs. [^|^ it was suggested to use a soft gluon exponentiation 
model to extend the perturbative calculation into this region of large minijet multiplicities. In 
the following we show that existing data on multijet events at the Tevatron allow one to 
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test and refine this model for multiple minijet emission. In particular one can experimentally 
determine "good" choices for the factorization and renormalization scales, information which can 
then be used to more reliably predict minijet emission in other processes, at the Tevatron or at 
the LHC. 

Consider dijet production in pp collisions at the Tevatron. The next-to-leading order (NLO) 
QCD corrections to the cross section for producing two or more jets, 0"2,inciE|, is available in the 
form of a full NLO Monte Carlo program, JETRAD [|l^]. The same program also provides the 



tree level cross section, a^, for three-jet production, within an arbitrary phase space region. In 
the following we use the CDF cone algorithm to define jets [|ri|], with a radius of Rq = 0.7 in 
pseudorapidity-azimuthal angle space; a cluster of partons of transverse energy Et > 20 GeV 
and pseudorapidity \ri\ < 4.2 is defined as a jet. (Variations of the Et threshold will be considered 
later.) Following the analysis of the CDF Collaboration we study two- and three-jet events 
with total jet transverse energy, 

Y,Et> 420 GeV , (1) 
and invariant mass of the multijet system, 

mjets > 600 GeV . (2) 



In addition the scattering angle |T2|, 9*, of the highest Et jet in the multijet center of mass frame 
must satisfy 

|cosr|<^. (3) 

Within these cuts (T2,inci depends weakly on the minimal jet transverse energy since the two-jet 
inclusive events are largely defined by the two hardest jets, which must have transverse energies 
Eti, Et2 ^Y.Et/2 > 210 GeV. The three-jet cross section. 



^ Here and in the following, the "cross section for n-jet inclusive events", cr„,inci = Z]fc>n c"fe jets, 
directly corresponds to the rate of such events. No jet multiplicity factor is included in its 
definition. 
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(T3{ET,min) = / ^^TSTTT^ , (4) 



on the other hand, is a steeply falhng function of the transverse energy threshold, Ej- j^i^. For 
sufficiently low threshold one will eventually reach a region with cr3(-Er,min) > cr2,inci; clearly, 
the interpretation of as the cross section for either three jet inclusive or three-jet exclusive 
events is not tenable in this region. The unphysical relation (T3 > o"2,inci is a sign that fixed order 
perturbation theory is breaking down and that multiple gluon emission needs to be resummed. 

For small i?T,min, soft gluon emission from the hard dijet production process will dominate, 
and, analogous to soft photon emission, one may show that this soft-gluon radiation approxi- 
mately exponentiates when the soft gluons go unobserved [0. Here we consider a phenomeno- 



logical model which assumes that the analogy to multiple soft photon emission can be taken 
further, namely, that the probability P„ for observing n soft jets beside the two hard jets of the 
basic hard scattering event is given by a Poisson distribution. 



with 



Pn{n) = ^e-\ (5) 



n = n(ET,mm) = — ^ / dET3 ■ (6) 

0"2,incl -^Et ,min GjH/'Y'^ 



We will call this model the "exponentiation model" in the following The exponentiation 

model has a number of appealing features: 



1. By construction the cross section for two jet inclusive events is given by cr2, 



incl- 



2. cr2,inci-Po gives the correct Sudakov suppressed rate for two jet exclusive events [0; this is 
not surprising since Pq is the probability for not seeing any additional jets, and here soft 
gluon exponentiation can be proved. 

3. For sufficiently large -ET.min, when <Js{ET,mm) ^ o"2,inci, the three jet rate, i.e. the cross 
section for events with one soft jet, is given by 

0"2,incl Pi ~ 0-2,incl fl = (y^{ET,rmn) , (7) 



and thus reproduces the perturbative result at 0{al). From a perturbative point of view 
o"2,inci-pL = o"3(-ET,min)(l + (^(^s)), thus, o"2,inci-Pi has the Same level of accuracy as the tree 
level calculation of cr3(i?T,min)- 

4. The full angular and transverse momentum information contained in the tree level result 
is retained in the estimated multijet emission probabilities. 

5. The probabilities Pn all remain finite at small i?T,min, with < P„ < 1. This renders the 
exponentiation model superior to the use of e.g. cts/ (cr2,inci — o"3) or crs/a2^inci as an estimate 
for the ratio of three-jet to two-jet exclusive or inclusive rates. 

6. Only the hard scattering cross section (here c'"2,inci) and the cross section for emission of 
one additional soft parton (here a^) need to be known perturbatively. Thus, the model can 
easily be applied to more complicated processes @J31- 

For the exponentiation model to be useful, it is necessary that minijet multiplicities in hard 
scattering events at least approximately follow a Poisson distribution. Recently the CDF collab- 
oration has published results on multijet production within the cuts of Eqs. (|T]-§). Out of a total 
of Ntot = 1874 events, A^o = 345 have exactly 2 jets with Et > 20 GeV, Ni = 612 have 3 such 
jets, and the number of events with n > 2 minijets, in addition to the hard dijet system, are 
N2 = 554, N3 = 250, = 88, N5 = 21, and A^e = 4 ||. For a Poisson distribution the average 
multiplicity of minijets in the CDF event sample, {rijgts ~ 2) = 1.57 ± 0.03 (where the error is 
statistical only), and the values for n{Nn/Ntot) extracted from the fraction of events with a fixed 
number of minijets, via the relation 

N 

fn = ^ = Pn{n) , (8) 

should all agree. The results of these different extractions of n are compared in Fig. |l]. A 
Poisson ansatz for the jet multiplicity distribution works surprisingly well. For larger minijet 
multiplicities {n > 4) the observed rates fall more and more below a Poisson distribution. This 
can be understood qualitatively, however, as an effect of the reduced phase space (in r], (p) which 
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FIG. 1. Poisson parameter n extracted from the fraction of events with n "minijets" in the 
CDF multijet sample according to Eq. (H). Errors are statistical only. A Poisson distribution 
of the minijet multiplicity would yield the dashed line at n = 1.57, the mean observed minijet 
multiplicity. No solution to Eq. exists for = 2 since the observed fraction, /2 = 0.296±0.011, 
is larger than the maximal Poisson probability for two events, P2{n = 2) = 0.271. The dashed 
line gives P2{n = 1.57) = 0.257 which is 3.5a below the measured f2- 

is available for additional minijets due to the finite cone size of jets. Also, one should note that 
a systematic error of 15% only on all Nn, (in addition to the statistical errors shown in Fig. |1|) 
would suffice to render a Poisson fit perfectly acceptable. 

Since the CDF data are well described by a Poisson distribution of minijet multiplicities we 
may now compare the measured average multiplicity, n = 1.57 for ET,mm = 20 GeV, with the 
perturbative QCD prediction given in Eq. (||). For this purpose we have calculated the total 
multijet rate, within the cuts of Eqs. (|l]-0) and 10% Gaussian smearing of jet momenta p to 



simulate detector effects. Using the JETRAD program |10] we find o"2.inci = 33 pb. Here we have 



used CTEQ 4HJ parton distribution functions [15| and the factorization and the renormalization 
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scale have been set to Hf = IJ^r = J2Et/4:, with as evaluated at two-loop order and as{Mz) = 

0. 116 as required by the evolution of the parton distribution functions. 

Since the 3-jet cross section cr3(i?T,min) is evaluated at tree level only, a variation of scales leads 
to substantial uncertainties here. We have analyzed four different choices, (i) fJ^R = ftp = J2 Et/4:, 
which is the same choice as for the 2-jet inclusive rate. Such a large factorization scale is not 
entirely physical, however. Collinear initial state radiation with transverse energies between 
ET,mm and the hard scale, I] -Et/4, is generated explicitly in terms of the third jet. When 
choosing a large factorization scale, this emission is considered twice, in terms of the third 
jet and via the evolution of the parton distribution functions. A factorization scale which is 
tied to the transverse energy of the third parton avoids such double counting. Similarly, a small 
renormalization scale may better match the small parton virtualities which appear in the emission 
of soft gluons. This motivates the following set of scales: (ii) fiR = fip = -E't.s, where Erp.s is the 
transverse energy of the softest of the three partons, and, finally, fip = ^Et,3 with an factor 
in the calculation of as which is given by 

al = llas{^ET,i) , (9) 

i=l 

1. e. the transverse energy of each final state parton is taken as the relevant scale for its production. 
Here we have used overall scale factors (iii) ^ = 1 and (iv) ^ = 1/2. 

Results for these four scale choices are shown in Fig. 0. For Ex^mm = 20 GeV, estimates for 
the mean minijet multiplicity n(i?T,min) = cr3(-E'T,min)/<72,inci vaij between 0.8 and 2.2 and thus 
bracket the CDF value of 1.57. From Fig. |^ it is clear, however, that a scale choice tied to the 
transverse energy of the soft jets is preferred by the data. The single CDF data point is not 
sufficient, of course, to optimize the scale choice: an analysis of the observed minijet multiplicity 
as a function of the minimal jet transverse energy would be needed. 

We conclude that the exponentiation model provides a good description of existing data 
on minijet multiplicities in hard scattering events at the Tevatron. For a further quantitative 
comparison, the ET,mm dependence of the multijet rates would be most useful. The freedom in 
choosing the renormalization and factorization scales for the tree level 3-jet cross section should 
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FIG. 2. Ratio of the tree level 3-jet cross section to the NLO cross section for 2-jet inclusive 
events within the CDF acceptance cuts [§. The cross section ratio n = cr3(i?T,min)/o"2,inci, with 
c"2,inci = 33 pb, is shown as a function of the transverse energy threshold, i^T.min, of the third 
jet. Results are given for four different scale choices in a^: fiR = ftp = J2Et/4: (solid line), 
I^R = I^F = Et^3 (dashed line), and fip = ^-Et.s, «s = IlLi «s(^-^r,i) with a scale factor ^ = 1 
(dash-dotted line) and ^ = 1/2 (dotted line). The CDF value for the average minijet multiplicity, 
n = 1.57, is given by the diamond. 

allow for a reasonable parameterization of the data, as a function of i?T,min and the parameters 
of the hard event, provided the approximate Poisson multiplicity distribution of jets remains 
valid beyond the phase space region probed in Ref. Finally we note that a comparison with 
tree level results for 4 or 5 jet cross sections [0, while complementary, will be subject to very 
large theoretical uncertainties. Just like in the exponentiation model, there is a strong scale 
dependence of the tree level n-jet cross sections (o"„jets ~ «"(/^k))- In addition one needs to use 
the tree level cross sections at the edge of the validity range of fixed order perturbative QCD 
and similar to the case of the 3-jet cross section discussed before, calculated n-jet rates will 



become unphysically large for sufficiently low £'r,min- The exponentiation model addresses this 
last problem and should, thus, provide a valuable alternative in analyzing multijet emission in 
hard scattering events. 
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